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SUMMARY 

A new {1, 2}-order theory is proposed for the linear elasto-static analysis of laminated 
composite plates. The basic assumptions are those concerning the distribution through the 
laminate thickness of the displacements, transverse shear strains and the transverse normal 
stress, with these quantities regarded as some weighted averages of their exact elasticity 
theory representations. The displacement expansions are linear for the inplane components and 
quadratic for the transverse component, whereas the transverse shear strains and transverse 
normal stress are respectively quadratic and cubic through the thickness. The main distin- 
guishing feature of the theory is that all strain and stress components are expressed in terms 
of the assumed displacements prior to the application of a variational principle. This is 
accomplished by an a priori least-square compatibility requirement for the transverse strains 
and by requiring exact stress boundary conditions at the top and bottom plate surfaces. 
Equations of equilibrium and associated Poisson boundary conditions are derived from the 
virtual work principle. It is shown that the theory is particularly suited for finite element 
discretization as it requires simple C°- and C '-continuous displacement interpolation fields. 
Analytic solutions for the problem of cylindrical bending are derived and compared with the 
exact elasticity solutions and those of our earlier {1, 2} -order theory based on the assumed 
displacements and transverse strains. 



INTRODUCTION 

Designing of aerospace and ground-vehicle structures with thick-section organic-matrix 
composites is necessarily associated with the analytical modeling of the structural response 
and the prediction of failure under service loads. For such applications, viable analytical 
models are those that can properly account for transverse shear and transverse normal 
deformations. These effects, which are negligibly small in thin laminates, can be significant 
when the laminate thickness and the wavelength of loading are the same order of magnitude. 
Moreover, the reduced stiffness and strength in the transverse shear and transverse normal 
material directions can contribute to significant deformations and matrix dominated failure 
modes such as delamination and transverse cracking. 


1 Present address: NASA Langley Research Center, Computational Mechanics Branch, 
Structural Mechanics Division, Mail Stop 240, Hampton, Virginia 23665-5225. 
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The majority of finite element methods for composite laminates are based upon the {1,0}- 
order theories (commonly referred to as lst-order theories). These account in some average 
sense for transverse shear but ignore the effects due to transverse normal deformations. The 
use of three-dimensional elements, including all three-dimensional effects, are generally 
computationally prohibitive even on a supercomputer, unless employed in some global-local 
fashion in a relatively small domain of interest. Although post-processing methods have been 
devised leading in some cases to adequate three-dimensional stress and strain recovery from 
the {1,0} -order theory computations, the static and dynamic response predictions obtained 
with this order of approximation are often inadequate for thick laminates (e.g., refs. (1-7) and 
references thereof). 

Numerous efforts to generate laminate theories of higher-order that also include transverse 
normal deformations failed to produce sufficiently accurate and computationally suitable 
formulations to model thick-section composites. Recently, Tessler (ref. 1) proposed a {1,2}- 
order theory for the elasto-static analysis of homogeneous orthotropic plates which was 
demonstrated to be ideally suited for finite element approximations. Subsequent developments 
included the analytic and finite element analyses of elastic beams (refs. 2-4), an extension to 
the general {1,2} -order orthotropic shell theory (ref. 5), and an extension to laminated 
composite plates for elasto-statics (ref. 6) and elasto-dynamics (ref 7). These theories assume 
linear inplane displacements and a quadratic transverse displacement across the thickness. 
The key aspect in the approximation is the independent expansions for transverse strains 
which allow exact stress boundary conditions at the top and bottom plate faces to be satisfied. 
These strains are also least-square compatible with those derived directly from the displace- 
ment assumptions. 

In this paper we propose a new { 1,2} -order theory for the linear elasto-static analysis of 
laminated composite plates following the methodology established in (ref. 6). The present 
formulation departs from (ref. 6) in that the transverse normal stress is independently 
expanded across the laminate thickness instead of the transverse normal strain. The transverse 
shear strains are taken to be parabolic satisfying traction-free boundary conditions. The 
transverse normal stress is assumed to be a cubic function through the thickness, also 
satisfying the consequence of zero shear tractions at the top and bottom laminate surfaces; 
specifically, the vanishing on those surfaces of the transverse normal stress gradient taken 
with respect to the thickness coordinate. These latter conditions are exact according to three- 
dimensional elasticity theory when the body and inertial forces are neglected. The resulting 
equilibrium equations are associated with exclusively Poisson boundary conditions. As its 
predecessor theory, the present theory requires only C°- and C'-continuous displacement 
interpolation fields and thus lends itself well to the development of simple, robust, and 
computationally efficient plate formulations similar to Mindlin-type elements (refs. (9-11)). 

In assessing the accuracy of the proposed theory, we resort to the problem of cylindrical 
bending of an infinite laminated plate for which an exact elasticity solution is available (ref. 
8). We compare the results of the present theory for several lamination patterns and span-to- 
thickness ratios with those of our previous {1,2} -order laminate plate theory (ref. 6) and the 
exact solutions. 
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membrane plate rigidities 
elastic stiffness coefficients for the k ,h ply 
membrane-bending coupling plate rigidities 
bending plate rigidities 
transverse shear plate rigidities 

continuous functions with discontinuous derivatives at element interfaces 
continuous functions that are discontinuous at element interfaces 
plate thickness 

inplane and transverse shear force resultants 
bending moment resultants 
applied transverse load at z=h 
midplane displacements along x and y axes 
Cartesian displacement components 
components of the transverse displacement 
inplane and transverse coordinates 
variational operator 
strain and curvature components 
bending rotations about x and y axes 
dimensionless thickness coordinate 
stress components 


{1, 2} -KINEMATIC PLATE ASSUMPTIONS 

In a laminated composite plate constructed of perfectly bonded orthotropic plies whose 
transverse constitutive properties do not differ appreciably, the displacement vector u = (u x> 
u y , u z ) can be approximated with functions that vary continuously across the total thickness. 
The lowest-order displacement expansions accounting for transverse shear and transverse 
normal deformations involve linear thickness variations for the inplane components u x and u y 
and a quadratic one for the transverse displacement u z (ref. 6) (i.e, {l,2}-order theory): 
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u x {x, y, z) - u(x, y) + hl,Q y (x, y ) , u y ( x, y, z) - v(x, y) + h t, Q x (x, y) 
u(x,y, z) - w(x,y) + Z,w l (x,y) + (^ 2 * * - 1/5) w 2 (x,y) 


(1) 


where ^=z/h€[-l, 1] is the dimensionless thickness coordinate and £=0 identifies the 
reference midplane position; w(x,y) and v(x,y) are the midplane displacements along the x and 
y axes, and 0,(x,y) and 0 y (x,y) are the rotations of the normal about the x and y axes (see fig. 
1). The Wj and w 2 variables can be interpreted as the normalized strain and curvature in the 
thickness direction: 
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The 3-D Hooke’s law for any kth orthotropic ply is expressed in the mixed form 2 
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(2.1) 


where C ij (k) are the elastic stiffness coefficients corresponding to the x-y coordinates (ref. 12). 


The components of strain and stress can be expressed in terms of the plate strain and 
curvature variables that are independent of the thickness coordinate: 


2 The variables superscribed with the k index are ply-dependent, whereas those without the 

k index are some average representations across the laminate thickness and thus are independent 

of the individual ply properties. 
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(3) 


e/ ” [e«. V “ K’ V + v ^ ] 

K / - K„> V K «’ K *>J - f 0 >^ e ^> e *, + e J 

y r -rv , V 1 - [w + 0 , w + 0 ] 

The inplane strains, obtained from eq. (1) in accordance with elasticity theory, are 
e „ - e x 0 + zK »’ e » - e ,o + z V Y^-Y „ + 


The present theory departs from that in ref. 6 only in the manner in which the transverse 
normal strain and stress are developed. Here, we independently expand the transverse shear 
strains and the transverse normal stress. The motivation for the latter assumption stems from a 
careful examination of a series of 3-D exact elasticity stress solutions (ref. 8) for a cylindrical 
bending problem (see fig. 2). Figure 3 shows 3-D solutions for the transverse shear and 
transverse normal stresses corresponding to the four distinct laminations: [0], [30/-30]s, 

[0/90]s and [0/90]. It is quite revealing that the form of the transverse normal stress remains 
relatively consistent from lamination to lamination whereas the transverse shear stress changes 
dramatically. These elasticity results clearly suggest that an independent, continuously varying 
expansion for the transverse normal stress should be an improvement over an analogous 
expansion for the transverse normal strain as used in (ref. 6). On the other hand, it is unlikely 
that one would gain any advantage by expanding the transverse shear stresses in a continuous 
manner across the laminate over that concerning the transverse shear strain expansions (ref. 6) 
since the solutions for these stresses (fig. 3) cannot be reproduced with a relatively low order 
approximation such as a parabola. 

The transverse normal stress and transverse shear strains are expanded independently across 
the total laminate thickness as 3 


- EoJ", Y,-, -Ey^" (i-x,y) (5) 

rt-o n-0 

The expansion coefficients G zn =a /n (x,y) and Y in =Y in (x,y) are determined by requiring the stress 
field to satisfy exact traction conditions at the top and bottom plate surfaces 

t f(x,y,±h) - o z Jx,y,±h) - 0 ( i-x,y ) 

along with the transverse strains to be least-square compatible with those obtained directly 
from elasticity theory using our approximations for the displacements (1): 


3 The present approach and that in ref. 6 are equivalent for homogeneous plates (ref. 1). 
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minimize J [ e'f - u zz \ 2 dz 

(7) 

minimize - (\, + M ii2 )] 2 dz O' - *,}0 


where the minimizations are performed with respect to the a zn and y jn variables. The resulting 
transverse normal stress and transverse shear strains can be expressed as 

- <U e „’ K o) + K „)^ - ^ 3 /3), Y 12 - j(l - ^ 2 )y, 0 (i-x,y) (8) 

Application of the virtual work principle results in the lO^-order plate equations of 
equilibrium and associated Poisson boundary conditions. The principle also serves as a 
variational framework for finite element approximations. The resulting two-dimensional 
variational statement has the form 

[N t 5e o + M t & k o 

where 5W, denotes the virtual work 
vectors of the plate stress resultants 
constitutive relations 


where A=[A jj ], H=[B y ], Z>=[D tj l and 

REMARKS ON FEM APPROXIMATIONS 

The variational principle (9) provides a convenient framework for developing efficient plate 
bending elements for the analysis of thin and thick composite laminates. In eq. (9), the 
kinematic variables u, v, w, 0 X and 0 y possesses spatial derivatives that do not exceed order 
one, thus requiring only C°-continuous finite element trial functions; w, and w 2 do not have 
spatial gradients thus needing only C 1 approximations. 

The simplest laminate plate element, which is also the most desirable from the standpoint 
of adaptive mesh refinement, is a three-node anisoparametric triangle (ref. 6). The element 
employs linear (3-node) parametric functions for u, v, 0, and 0 y , an anisoparametric quadratic 
field for w and uniform fields for w, and w 2 . The latter variables only contribute two degrees- 
of-freedom per element; in static problems, they can be conveniently condensed out at the 
element level using static condensation. 


+ Q T Sy t ]dxdy-W'(u,v,w,Q x ,Q y ,w v w 2 ) - 0 (9) 

of external forces; A^={ N fj } , A/={Mjj} and 0={Qi} are 
which are related to the plate strains (eqs. (3)) via the 
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G=[Gj|] are the plate constitutive matrices. 
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RESULTS AND DISCUSSION 


The present theory will be qualitatively assessed with respect to its application to the 
problem of cylindrical bending of an infinite carbon/epoxy laminate subjected to a sinusoidal 
transverse pressure q + = q 0 sin( 7 tx/L). The exact elasticity solution for this problem was first 
determined by Pagano (ref. 8). We shall consider an orthotropic [0] laminate, two symmetric 
laminates — a cross-ply [0/90], laminate and an angle-ply [30/-30] s laminate — and an 
antisymmetric cross-ply laminate, [0/90]. The ply material properties are taken as 

E l - 25x1 0 6 psi, E t - 10 6 psi, G LT - 0.5x10 s 

G n - 0. 2x\0 6 psi, v LT - V-jt m 0.25 

where L and T denote the longitudinal and transverse ply material directions, respectively. 

As noted previously, this theory differs from ref. 6 only in the manner in which the 
transverse normal stress and strain are approximated. Thus, it is only natural to expect that the 
differences in the predictions with the two theories will mostly affect these transverse normal 
variables. Indeed, the displacement predictions by both theories are virtually identical, and so 
are the inplane and transverse shear stresses and strains. 

Figure 4 depicts the percent error in the maximum deflection computed at the midplane of 
the laminate versus the length-to-thickness ratio L/2h. The results, which pertain to the 
present and ref. 6 theories, clearly show that as far as the deflection predictions are con- 
cerned, the engineering accuracy (i.e., error < 5% ) is attained for laminates with the ratio 
L/2h > 4. These results demonstrate the adequacy of the overall plate stiffness representation. 

Figures 5 through 8 compare the stress and strain thickness distributions for thin (L/2h=40) 
and thick (L/2h=4) laminates obtained with the present theory (designated as HOT-S), our 
previous {1,2} -theory (HOT-E, ref. 6), and the exact solutions obtained in this effort using the 
3-D elasticity approach (ref. 8). Figure 5 shows the transverse normal stress distributions 
across the laminate thickness. Note that while both HOT-E and HOT-S are nearly equally 
accurate in the [30/-30] s lamination predictions — with HOT-E exhibiting a slight discontinu- 
ity at the interface between 30 deg. and -30 deg. plies — the HOT-S prediction for the cross- 
ply [0/90] laminate is much superior. The transverse normal strain distributions are depicted 
in fig. 6, where the present theory exhibits proper discontinuous character at the ply interfac- 
es. Except on the bottom surface, both HOT-E and HOT-S produce appreciable errors in the 
[0/90] thick laminate. 


The transverse shear stress (t xz ) distributions are compared in fig. 7. The HOT-S and HOT- 
E stresses are computed by integrating appropriate inplane stress gradients in the 3-D 
equations of equilibrium (ref. 6). Throughout, these predictions compare well with the exact 
elasticity solutions. Figure 8 depicts thickness distributions of the normal stress (o xx ). These 
show that for the [30/-30]s thick laminate the stresses in the outer plies are underestimated. 


1423 



This is naturally the result of linear assumptions for the inplane displacements (1). Interesting- 
ly, the results for the [0/90] thick laminate are excellent. 

CONCLUDING SUMMARY 

In this paper we have discussed a new {1, 2} -order theory for the elasto-static analysis of 
laminated composite plates in which the assumed variables are the components of the 
displacement vector, the transverse shear strains and the transverse normal stress. The 
transverse stresses satisfy exact stress boundary conditions at the top and bottom plate 
surfaces. The virtual work principle produces a set of equilibrium equations and associated 
Poisson boundary conditions. As was the previous theory, this theory is particularly suited for 
finite element approximation with simple C°- and C '-continuous displacement interpolation 
fields. The analytic solutions - obtained for a wide range of laminations and thicknesses - 
showed that the theory is applicable to thin and thick laminated composites and has some 
advantages over our previous {1, 2} -order theory. 
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